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Abstract. This article reports an investigation of the effect oi soion: rammrpmpte^iire on 
the motion of an artificial satellite. The theory has been applied to the orbit of the Vanguard 
I satellite, and is found to produce significant perturbations in the perigee height of that 
satellite. In the case of a satellite with a large ratio of area to mass the major terms introduced 
by solar radiation pressure can reduce the perigee distance at the rate of 1 to 2 km per day, 
so that the lifetime of the satellite will become considerably shorter than it would be without 


this effect. 

Introduction. The perturbing effects of solar 
radiation pressure on satellite orbits have been 
generally considered by celestial mechanicians to 
be negligible. However, the difference between 
observed and theoretical values of perigee height 
for the Vanguard I satellite have suggested a re- 
examination of radiation pressure as a possible 
source of the discrepancy. We have carried out an 
investigation of the radiation-pressure perturba- 
tion and developed an analytical theory of 
resonance for the case in which the perigee 
follows the motion of the sun. It is found that 
when resonance conditions are nearly satisfied 
the radiation pressure may produce substantial 
orbit perturbations over an interval of several 
months. For a satellite with a large ratio of 
area to mass the major terms introduced by 
solar radiation pressure can reduce the perigee 
distance at the rate of 1 to 2 km per day, so 
that the lifetime of the satellite becomes con- 
siderably shorter than it would be without this 
effect. Radiation pressure is also found to produce 
significant perturbations in the orbit of the 
Vanguard I satellite which remove the principal 
discrepancy between theory and observation in 
the analysis of the Vanguard I orbit [ Musen , 
Bailie, and Bryant, 19G0]. 

The terms of long period make the principal 
contribution to the orbit perturbations. In the 
present investigation we have therefore neglected 
the effect of the earth’s shadow which has the 
period of the satellite’s mean motion. The 
shadow effect does not change the nature of 


the perturbations, although it produces a change 
in their amplitudes. We have also neglected all 
other periodic terms having the mean anomaly 
of the satellite in the argument. 

The vectorial method is used in the develop- 
ment to obtain several kinematical relations. 
The equations for the scalar elements are deduced 
from the equations for the vectorial elements. 
All elements except the semimajor axis contain 
long-period terms, but the semimajor axis is 
affected by short-period variations only and is 
therefore not subject to substantial perturba- 
tions. 

Perturbations in the orbit plane. Let us 
designate the gravitational constant by G , the 
mass of the earth by M and its equatorial radius 
by p. Let R be the unit vector directed along 
the normal to the orbit plane, P be the unit 
vector directed from the center of the earth to 
perigee, and let Q = R X P. The position and 
the velocity vectors of the satellite are repre- 
sented by r and v, respectively; the radius vector 
of the satellite will be designated by r, the true 
anomaly by /, the mean anomaly by l\ and the 
other elliptic elements will be designated, using 
the standard notations, by co, i, e, a, n. The 
mean longitude of the sun on the ecliptic will 
be designated by X', and the mean motion of the 
sun will be designated by n'. Let c be the inclina- 
tion of the equator to the ecliptic, i the unit 
vector directed from the center of the earth 
toward the vernal equinox, k the unit vector 
normal to the earth’s equator, and let j = k X i. 


1391 


V 


1392 


PETER MUSEN 


Finally, let u° be the unit vector directed from 
the center of the earth toward the sun. Neglecting 
the eccentricity of the earth’s orbit, wc can put 

u° = i cos X' + j cos e sin X' 

+ k sin € sin X' (1) 

The solar radiation pressure will be 


F = Fu 


F < 0 


( 2 ) 


F = constant for a spherical satellite; for the 
satellite of nonspherical form we must make a cer- 
. tain assumption about the average value of F and 
use the average value in the development. We 
shall deal with the vectorial element 


g = cP 


(3).... 


and with the time variation of this vector with 
respect to a system of coordinates rigidly con- 
nected to the osculating orbit plane. In other 
words, we shall deal with the motion of P in 
this plane. For that purpose it is convenient 
to use Herrick’s equation [ Herrick , 1948] 

GM dg/dt ~ T-¥ (4) 


where 


2rv — vr — r-v/ (5) 


r and v in this case are given, not with respect 
to our inertial system, but with respect to a 
system rigidly connected to the osculating orbit 
plane; The notations rv and vr represent the 
dyadic products; in other words, they are the 
products 

column vector- row vector 

and I is the planar idemf actor (the planar unit 
matrix). Substituting 


-P 


r = Pr cos / + Qr sin / 
VGM sin / 


V o(l — 

+ Q 


2\ 

e ) 

V GM (cos / + e ) 
y/ a(l — e) 


into consideration, wc deduce 

r - 

+ 5 sin 2/)PP 

+ sin 2/QQ 
a a 

+ | - (-H- e cos / \ 

l_a 

+ i cos 2/) + 1 - c 2 ]pQ 

" [i a ( ‘ “ "■ 2/) 

+ 1 - e 2 ]Qp| 

In order to separate the long periodic part in 
T from the short periodic one we must develop 
the coefficients of the dyadic products into 
Fourier series with respect to the mean anomaly 
and retain only the constant terms in this 
development. Using Cayley’s tables, we find 
that the constant parts of the coefficients of 
PQ and of QP are- equal to ifc|(l — c 2 )- The 
constant parts of the other coefficients are equal 
to zero. Thus, we have for the long, periodic 
part of T 


I = PP + QQ 
into (5), and, taking 

GM = na 


( 10 ) 


[r] 




Vl - e 2 (PQ •- QP) 


but 


R X I = R X (PP + QQ) = QP - PQ 
and we have a simple relation 

[r] = - 


2 Vl - e R X / (11) 


( 6 ) 

(7) 

( 8 ) 

(9) 


Substituting (3) and (11) into (4) we deduce 
for the long periodic part in eP 


d(eP) 

dt 


3 na 
~2 


Vi - 


GM 


R X F 


( 12 ) 


We see that the time variation of eP is normal 
to the component of F lying in the orbit plane. 
The vector P in the moving system of coordinates 
is affected only by the rotation about R. The 
angular velocity of rotation 

R dw/dt 



Thus 
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os 


3 F 


'na y/ 1 — e 


It follows from (12) and (13) 


' de 


I 5 p + e 


dir 

dt 


3 ha y / 1 — 
2 


GM 


R X F 


de _ 3 nd_ 

dt “ 


Vi~- 


(?M 


Q F (14) 


dt 2GM 

{4 cos 2 (i/2) sin 2 ( e/2 ) cos (co + XX + X') 

+ cos 2 (i/2) cos 2 (e/2) cos (co + XX ~ X') 

4s in 2 (t*/2) cos 2 (e/2) cos (co — XX 4 X') 

+sin 2 (i/2) sin 2 (e/2) cos (co — XX — X') 

— \ sin i sin e cos (co + X') 

+ 2 sin i sin e cos (co — X')} (10) 


dir _ 


3 na 2 \/T 
2~ 


GM 


PF 


dt 

Substituting 

P = i[+cos 2 (i/2) cos (co 4 ft) 

4* sin 2 (i/2) cos (co — XX)] 

+ j[4cos 2 (i/2) sin (co + XX) 

— sin 2 (i/2) sin. (co — XX)] . 
4 k sin i sin co 

into (15) and 

Q = i[ — cos 2 (i/2) sin (co + XX) 

.. — sin 2 (i/2) sin (co — XX)] 

4 j[+cos 2 (i/2) cos (co + SI) 

— sin 2 (i/2) cos (co — XX)] 

+ k sin i cos co 


Perturbations in the position oj the orbit plane . 
(15) ” The angular velocity of rotation of the osculating 
orbit plane, considered as a rigid body, 


^T= 


na 


r RF 


Vi - c 2 a GM 
The part of r independent of l is 


(20) 


(16) 


and, consequently, the long periodic part of 
i is 


[*] = ~ 


nae 


vT 


PR - F 

GM 


( 21 ) 


The variation of R is caused only by the rotation 
of the osculating orbit plane, and we have then 
for the long periodic part in R 

dR/dt = -R X [ifcj 


(17) 


or 


into (14), and taking (1)— (2) into account, we 
obtain 


dR 

dt 


+- 
1 n 


na/eQ R - F 


2 GMVl - c 2 


( 22 ) 


de 

dt 


'VT^ 


and we conclude that the long-period time 
variation of R, caused by the radiation pressure, 
consists of the rotation of R about P with the 
{ 4cos 2 (i/2) sin 2 (e/2) sin (co 4 XX 4 ^0 angular velocity proportional to the cosine of the 


3 F 
2GM 


4 cos 2 (i/2) cos 2 (e/2) sin(co 4 XX — X') 
4s in 2 (i/2) cos 2 (e/2) sin (co — XX 4 X') 

4sin 2 (i/2) sin 2 (e/2) sin (co — XX — X') 

— \ sin i sin e sin (co 4 X') 

4 i sin i sin e sin (co — X')} (18) 


angle between the direction to the sun and R. 

This time variation is zero at the moment 
when R is normal to u°, in particular when R is 
normal to the ecliptic. From equation (22) we 
can also deduce the equations for i and XX. The 
angular velocity [t|j] can be decomposed into 
the geometrical sum of (1) the angular velocity 
di/dt of the rotation of the orbit plane around the 
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line of nodes and (2) the angular velocity, sin iw 
dQ>/dt of the rotation around the vector P sin co 
+ Q cos to. We have for the long-period variar 
tion in t and Q>, taking (21) into account, 


(P cos to — Q sin co) 


di 

dl 


+ (P sin co + Q cos co) -sin t 


. d£l 
dt 


nae 


'2 vT 


PR F 

GM 


We deduce from this last equation, after the 
scalar multiplication by P cos co — Q sin co and 
by P sin co + Q cos co and taking (P cos co 
Q sin co)*k = 0 and P k = sin i sip co into 
consideration, 

di 3 mx 2 cR-F cos co 

dt ~ 2 GMVi-7 

.dQ>_ 3 nacR' F sin co 

Sin 1 hi ~ 2 gm Vi - ? 

Taking 

R = i sin i sin — j sin i cos Q> + k cos i (25) 

and (l)-(2) into consideration, we have after 
some easy trigonometrical transformations 


+ sin i cos 2 (c/2) cos (co + & — X') 

— sin i cos 2 (e/2) cos (co — + X r ) 

+ cos i sin c cos (co + X') 

— cos i sin c cos (co — X') } (27) 

These equations could also be obtained in the 
standard way. 

The integration problem. If we assume that 
we do not have sharp resonance conditions, the 
first-order perturbations can easily be obtained 
by the integration of equations (18), (19), (26), 
and (27), assuming that the elements in the right- 
hand sides are replaced by their mean values. The 
oblateness of the earth produces the secular 
motions of the node and of the argument of the 
perigee, for which we have 


* (l ” to) 



Gk 2 p 2 ny[ — 

5.2. 
7 sin i 
A 

(23) 

to = too + 2/., 

a (1 *— 

2\2 
G ' 


3/c 2 p 2 7l COS 1 


(24) 

— £<!U 2/-. 2\ 2 

a (1 — e ) 

(t t ( 


The resonance case. The resonance case 
deserves special attention and a special treat- 
ment. The most interesting resonance occurs 
when the perigee of the satellite closely follows 
the sun, i.e., when the critical argument, in the 
terminology of celestial mechanics, is 


di' _ 3 Fnae 

dt ~ 4 GMVl - ? 

« { +sin i sin 2 (c/2) sin (to + + X') 

— sin i sin 2 (c/2) sin (to — SI — X') 

+ sin t cos 2 (c/2) sin (to + & — X') 

— sin t cos 2 (c/2) sin (to — + X') 

+ cos t sin c sin (co + X') 

.— cos i sin e sin (co *— X')l 

. . dQ> , 3 Fnac ’’ - 

sin i — = -tt 7 ~~ 2 

dt 4 0 M/V 1 ~ e 

• { +sin i sin 2 (c/2) cos (co + £1 + X') 

— sin i sin 2 (e/2) cos (co — Q> — X') 


co -i- ft - X' 

The term with this argument is the most in- 
fluential one in the development of the radiation 
perturbations in the case of Vanguard I. 

In order to simplify the exposition we shall 
adopt the system of units in use at the Vanguard 
Computer Center in the computation of the 
general oblateness perturbations. We put 

G = 1 M — 1 p - 1 


and we use the system of coordinates rotating 
uniformly in the equatorial plane with angular 
^ velocity equal to the mean motion of the sun in 
the ecliptic. In our exposition wc follow the 
line of thought laid down by Brown in his 
planetary theory [Brown and Shook, 1933], 
modifying the form of the disturbing function 
to serve our purpose, and, as Brown did in the 
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planetary case, we will reduce our problem to 
an equation similar in form to the equation of 
the motion of a mathematical pendulum. The 
modified canonical set of Poiucard 


Xl = (1 - V J — e") 

?/i = - w — ft + A' 
x a = y/a V 2 — l + «> + ft — A' 
xa 


(28) 


= VoO - ?) (1 - COS i) 

t/ 3 = — ^ + x / 


is especially convenient here. Phe corresponding 
Hamiltonian function consists of several parts: 


1. The part 


2 a 


is contributed by the elliptic motion of the 
satellite. 

2. The part 


n f \/a( 1 — e) cos i 
is generated by the motion of the sun. 

3. The secular part 

k 2 (l — I sin 2 i) 

3/, 2x3/2 

CL \\ C ) 

is produced by the oblateness of the earth. 

4. The periodic part 

— § Fae cos 2 (i/2) cos 2 (*/ 2) 

• cos (co + & " 


+ 


K* ■ - f Foe cos 2 (t‘/2) 


3/1 J^ 3/2 

a (1 — e J 
• cos 2 (e/2) cos (co + & — X') 
or, taking (28) into consideration, 

1 


i r i _ 3x 3 

+ x 2 3 L(x 2 - X,) 3 (*» - xi) 4 

Xl) R l 


+ k-Ji cos 7/, 
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(30) 


2(x 2 ~ 
where we put 

k 2 K — ~- '£F(\c cos 2 (i/2) cos 2 (*/2) > 0 
We have 


A') 


comes from the development of the disturbing 
potential — F r into a trigonometric series. 

Thus, we have 

// = ~h n f \/fl( 1 6 ) cos t 

2a 




2x 2 


+ n'(x 2 — Xi — x 3 ) 


dx\ _ . dH 
dt 'dy i 


dy t 

dt 


dH 

dxi 


(31) 


The arguments yz and y z are not present in (30), 
and from 


dx 


dll 


_2 = « 0 
dt ■ dy 2 


dx 3 __ , dH _ ^ 
dJ "oifo 


it follows that 

x 2 = a 2 = constant 
x 3 = a3 = constant 
and we also have the integral of Jacobi 
II = h 

If the satellite is approximately of the size and 
the form of Vanguard I, the perturbations in 
the eccentricity are small, even if the resonance is 
a sharp one. The ratio Fjkz in this case is also 
small, and the coefficient K might be considered 
to be invariable. Let us introduce a new variable 
£ instead of Xi, and a new independent variable 
T instead of t } by putting 

Xi = ai + £ (32) 

Vh dt « dT (33) 

where ai is a constant to be determined. 

Equations 31, if transformed to the new 
variables, becomes 


di 

dT 


+ 


d W 
dy i 


dyx 

dT 


dW 
' d£ 


(34) 


(29) 


where 
W 


" (‘- e*. 


+ k a g + K cos Vi 

OCX 1 


(35) 
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. r r i jga 

• a* L(a z , «1 Y (« 2 — «i) 4 

‘ + 3W _1 

■■■■ / 2(a 2 -a 1 ) 5 J 

Let us determine the constant « t in such a way 
that the condition 


v = k 2 K<t>a\ * 


(42) 


(36) 


and 0 < < 1 is the constant of integration. 

We deduce 


dy } 

dT 


= — 2^ cn(vT , fc) 


(43) 


n' fc 2 |A = 0 

a«i 


(37) 


, is satisfied. This last equation says that in the 
first approximation the argument y i does not 
have a secular motion. In other words, to assert 
(27) is the same as to assert the existence of the 
resonant conditions. Equation^ (35) becomes 

: w 

and 


\k 2 *—5 £ + K cos y x (350 


dai 


- .. " da* > ° 

for moderate inclinations. We deduce from (34) 
and (350 


and, taking (32), (33), (42), and (43) into con- 
sideration, we obtain 

I K. 

Xi = ai + 2/cA/777-“cn(^7 T , k) 

>0n, a 

If the motion of yj is progressive, i.e., if yi has 
secular motion, then we deduce 

y x — tt — 2am(pT , k) 

= 7 r — 2am^~ V L ^ 

♦ 

and 

= W„,>" '' 


dT 

d]h 

dT 


— — K sin y x 






From tli esc two last equations we obtain 


dj/> 

dT 2 


K sin y, 


(38) 

(39) 

(40) 


*, * a. + | fc) 

And, finally, in the case of the asymptotic 
approach, 




» ?/ 1 - p T 

tg - = Cc 


k 2 K<j>a\* 


— > 0, £ — > 0, 


«1 


The integration of this equation is given in 
treatises on analytical dynamics, and it is not 
necessary to reproduce dt here. The following 
three cases exist 

1. The critical argument yi oscillates around 
180°; we have the case of libration. 

2. The critical argument has a secular motion. 

3. The critical argument approaches zero 
asymptotically. 

Despite the fact that the libration of the critical 
argument may be large, the perturbations in 
xi are always small. 

, In the case of libration wc have 

cos \yx = k sn (v7\ Jc) 

= ksn t, *) (41) 


• - Vi 

if r— » + 
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